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Abstract 

In this paper we consider an interacting particle system modeled as a 
system of N stochastic differential equations driven by Brownian motions 
with a drift term including a confining potential acting on each particle, and 
an interaction potential modeling the interaction among all the particles of the 
system. The limiting behavior as the size N grows to infinity is achieved as a 
law of large numbers for the empirical process associated with the interacting 
particle system. 



Introduction 

We consider a system of N(€ N — {0}) particles. From the Lagrangian point 
of view, the system is described by N random variable, X^(t) G M. d , t > 0, k = 
l,...,N, so that {X]fr(t),t £ R+} is a stochastic process in the state space 
(R d ,£> R d), on a common probability space (f2,.F, P). may describe the 
state of the k-th particle, e.g. its position. We consider the case of contiuous 
time evolution, i.e. the time evolution is described by a system of stochastic 
differential equations (EDSs) with additive noise 



dX k N (t) = f(X k N (t))+F N [X N (t)]{X k N (t) 



dt 

+ a N dW k (t), k = l,...,N. (1) 

In equation (pQ) the process {W k ,k = 1,...} is a family of independent 

standard Wiener processes, / : R d H ► R, and the functional Pjy is defined 

on A4-p(R d ), the space of all probability measures on R rf , and depends on the 
empirical measure 

1 - 

X N( t ) = xY, e x%(t)£M P (R d ). (2) 

k=l 
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By the empirical measure ([2]) we describe the system by an Eulerian ap- 
proach: the collective behavior of the discrete (in the number of particles) 
system, may be given in terms of the spatial distribution of particles at time 
t. 

Correspondingly, the measure valued process 

1 N 

X N : t G R + - X N {t) = jjY1 e x%(t) ( 3 ) 

k=l 

is called the empirical process of the system for a population size N. 
The trajectories are random elements of C([0, T], Ai-p(R d )), so that the 
distributions C(Xj^) of those processes can be considered as elements of 
M v (C([0,T],M<p(R d ))). 

Equation ([I]) might describe a system of N individuals whose movement 
is due to a stochastic individual component coupled with an interaction term, 
and an (individual) advection term. 

The individual dynamics The advection term / : ~R d -\ ► E may de- 
scribe the individual dynamics of a particle, which may depend on external 
information. Indeed we consider the following form for / 

f{x) = - 7 iW(x), (4) 

where 71 € R+, and the potential U : M. d — > M + G C 2 (R. d ) is a non negative 
smooth even function. From the modelling point of you the transport term 
([4]) mean to be "confining" potential: there are some external information 
coming from the environment which attracts the particle along the flow due 
to U. 

The interaction dynamics Fn it depends on the relative location of 
the specific individual X^(t) with respect to the other individuals, via on the 
empirical measure of the whole system. The interaction we consider is due to 
different phenomena: aggregation and repulsion. These two different forces 
compete but act at different scales. 

Aggregation act at macroscale and is modelled by a McKean-Vlasov in- 
teraction kernel 

G a :R d — ► M+. 

The interaction of the particle located in X^(t) at time t with the others 
is described by a "generalized" gradient operator as discussed in [U [14] acting 
on the empirical measure 

(VG a *X N (t))(xUt)). (5) 

Repulsion acts at mesoscale; the mesoscale is introduced as in [114 [14} [TB] 
via a rescaling of a referring kernel V\ 

V N (z) = N' 3 V 1 (N^ d z), 0E(O,1) (6) 
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The repelling force exerted on the k-th (out of N) single particle located 
at X N (t) dis given by 

N 

-^N^VV^N^fx^-X^t))) = -(VV N *X N (t))(X k N (W) 

i=l 

From ([7]) it is clear how the choice of (3 determines the range and the 
strength of the influence of neighboring particles; indeed, any particle interacts 
(repelling) with O (TV 1 "' 3 ) other particles in a small volume O (N~p). 

From ([5|) and (|7|). the advection interaction term F[Xn] is given by 

F[X N ](x) = l2 (V(G a -V N )*X N (t))(x), (8) 

with 72 G M + . 



The stochasticity The stochastic component in equation ([T|) may de- 
scribe both the lack of information we have about the environment or the 
particle itself and the need of each particle to interact with the others, so that 
they move randomly with a mean free path on (depending on N) unless they 
meet other particles and interact. 

By (pQ), ®) and © the system we study is the following 



dX N (t) = 7l VU(X N (t))+ l2 V((V N -G a )*X N )(X N (t)) 



dt 



+a N dW k (t), k = l,...,N. (9) 

In the case 71 = 0, the advection is due only to the interaction and the system 
become the following 



dX k N {t) = (VG a * X N (t))(X k N (t)) - (VV N * X N (t))(X k N (t))dt 

+a N dW k (t), k = l,...,N. (10) 

In previous papers [U [TTj the authors has focused their attention on 
the time evolution of the system (|10p . In particular they have analyzed the 
convergence of the system as the number of particles N increases to infinity. 
In [T3] a " law of large numbers" is presented while in [1] the authors have 
studied the existence and uniqueness of the solution of the PDE describing 
the time evolution of the limit system. 

In this work we focus our attention on the system ([9|) and extend to this 
case the results obtained in [Til E] • 



Notations and Hypotheses 

For some topological space S we denote by C™{S, M. d ) the space of m-times 
differentiable M d -valued functions on S with continuous bounded deriva- 
tives; C fe m (5,M) is abbreviated with C^{S) and C$(S,R d ) with C b {S,R d ). 
C™(S, M. d ) is equipped with the supremum norm. On R d x R d , (•,•) denotes 
the usual scalar product. 
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Ai-p(S) is the space of probability measures on S. This space is equipped 
with the usual weak topology. On the space A4-p(M. d ) the weak topology is 
generated by the complete metric 

H/x -v\\i= sup ((fij) - (v,f)), 
feHi 



where 



<A*,/>= / f(x)Kdx) feC b (R d ), 



and 



.(R d ) : sup < 1, sup lf{x) f l y)l < 1 



The metric di{n,v) := ||^ — u\\\ is also well known as bounded Lipschitz 
metric. 

For any 5- valued random variable Y we denote by C(Y) G M-p(S) its 
distribution. 

For some T G (0, oo), C([0,T], Mp(R d )) is the space of all continuous 
functions / = f{t), < t < T from [0, T] to M v (R d ), equipped with the 
metric 

p(f,g) = sup ||/(t)-s(*)lli- 

0<t<T 

For / G L 2 (R d ) we denote by 

/(A) = lim f — ) ' f e iXx f(x)dx 

a->oo\2irJ J{\x\<a} 

its Fourier transform. 

In connection with Fourier transforms we shall use the relations 



/ 



f{x)g(x)dx = f(\)~g(\)d\ /, ff e^(r), (11) 
JR d 

jr g (X) = (2nf 2 f(\)~g(\) f,g G L 2 (R d ), (12) 
Vf(X) = -iXf(X) feWi(R d ); (13) 

where 

= / (l + |A| 2 )|/(A)| 2 dA = ||/||i + ||V/||2<oo}. 

JR d 

Positive constants throughout the thesis are denoted by c\ , C2 , . . . ; if a constant 
depends on a quantity k, we denote it with c(k). 
By defining 

g N (x,t) = {V N *X N (t)){x), (14) 
Vg N (x,t) = {VV N *X N (t)){x), (15) 
h N (x,t) = (W N *X N (t))(x), (16) 

where 

W N (x) = x d N W 1 ( XN x) (17) 
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and W\ is a symmetric probability density denned on R d , that is 

W 1 (x) = W 1 (-x), (18) 

equation ([9|) becomes 

dX k N (t) = - 7l VU(X k N (t))+ 72 [(VG a *X N (t))(X N (t))-Vg N (X k N (t),t) 

+ a N dW k (t), k = l,...,N. 
We consider the further assumptions: 



dt 



V 1 (x) = (W 1 *W 1 )(x) = / W 1 (x-y)W 1 (y)dy 



and 

W 1 €W%(R d ) i.e. / (1 + |A| 2 )|VFi(A)| 2 (iA < oo. 
About the initial condition we suppose that 



sup E 

7V€N 



\x\X N (0)(dx) 



< oo, 



sup E 

tvgn 



\h N (x,0)\ 2 dx 



supE[||^(-,0)|||] <oo, 
NeN 



(19) 
(20) 

(21) 
(22) 



V N > 1, Xjv(0) = (X 1 {0),...,X N {0)) is independent of W k , k = 1, . . . , N 

(23) 

The assumptions about the interaction potential are the following: 

G a eC 2 (R d ,R + ) (24) 
is a symmetric function on R d supposed to be independent on N, 
Vn is supposed to be of the form V/v(x) = Xn^(xnx)-, where 

Vi G C 6 2 (M d ,IR + ) is a symmetric probability density on R d , (25) 



XN = N^ d , /3e(0,l). 



(26) 



It is clear that 



lim Vn = 5 , 

JV— >+oo 



where 5q is the Dirac delta function. 

VJV>1, E{|Xjv(0)| 2 } < +oo. (27) 
Let the confining potential be such that 

U £ Cl(R d ,R+)nC 2 (R d ,R+). (28) 

Then we consider the following possible assumptions on the parameter (3: 



a) 



e 0, 



, lim o";v = or 
d + 2 J 1 N^oo 



d 



d + 2 



,lV lim a N N^ d+2 y d ~ 1 = 0, 



b) 



/3 € 0, — — , lim a N = a 00 > 0. 

d + 2 / TV^oo 



(29) 



(30) 



A Law of Large Numbers 

In this section we derive a law of large numbers for the measure valued process 
Xn defined by ([2]) and ([3]), in the case of boundedness properties of the 
confining gradient U. In particular by following the approaches proposed in 
|12j and [15], we prove the existence of the limit measure for the sequence 
{£(Xjv)}7veN of distributions of {^OvIiVsN- 

We consider both the unviscous case, that is lim^v— »oo on = 0, and the 
viscous case, i.e. lhmv^oo = > 0. 

The procedure may be divided into three steps: 

i) relative compactness of the sequence £(Xn), N G N, which corresponds 
to an existence result of the limit C(X); 

ii) regularity of the possible limits: we show that the possible limits 
{X(t),t £ [0, T]} are absolutely continuous with respect to the Lebesgue 
measure for almost all t S [0,T] P — a.s.; 

iii) identification of the dynamics of the limit process: all possible limits are 
shown to be solution of a certain deterministic equation that we assume 
to have a unique solution. 

In the case limTv^oo <jn = we guess the limit dynamics and show that it is 
the weak limit of {X N (t),t £ [0,T]}. 



Relative Compactness 

The first step toward proving a law of large number for a measure-valued 
process is to obtain a relative compactness result for the sequence of em- 
pirical measure's distribution laws {£(Xjv)}jv e N associated to the system of 
stochastic differential equations. 

Theorem 0.1. // either ( fgg|) or fgQj) holds, under conditions p4\)-[2E\). the 

sequence {C(Xn)}n(z^ of distributions of the processes {XN(t),0 < t < T} 
associated to the system of stochastic differential equations (0) is relatively 
compact in the space M v (C([0,T],Mv(R d ))). 

We consider first some preliminary results regarding the martingale prop- 
erties of some processes. Up to now we suppose that all the hypotheses of 
Theorem 10.11 are satisfied. We remark that all the results are valid also in the 
case 7i = 0. 



For the seek of simplifying the notations, in the following calculations we 
set 71 = 72 = 1. 
Let 

A N (t) = [\x N (s),2(\Vg N (;u)\ 2 - Vg N (;u)(-VU(-) + (VG a *X N (u))(-)) 
Jo 

+ | - W(-) + (VG a * X N (u))(-)\ 2 )) + a 2 N \\Vh N (;u)\\ 2 2 du. (31) 
Lemma 0.1. The process 

M N (t) = \\h N (-,t)\\ 2 2 + A N (t)- [\x N (u),2\-VU(-) + (VG a *X N (u))(-)\ 2 )d u 

Jo 

- c l a 2 N tN^ d+2 ^ d - 1 
is a martingale. 
Proof. 

Because of (J23J), (HHJ, © and ©, by applying Ito's formula to 

1 ^ 

||M-,t)||! = <*iv(*),<?iv(-,*)> = T72 E ~X l N (t)), 



one obtains 

E[||M^)llll^] = \\h N (;s 
+ E 



2 

N 



f £ (^jvN, VG (- - • V 5 jv(-, u))d« 

m=l 



2 



s m=l 
JV 



E ^( u )' ( VV ^(- " X n("))) ■ Vg N (; u))d 

m=l 

2 r t N 



ii 



fc,m=l, k^trn 



\\h N (;s)\\l 



2 / (Xjv(ii),|V5iv(-,«)| 2 )d« 



E 



2 / (X A r(n),V 5J v(-,n)-(-VC/(-) + (VG a *X JV ( U ))(-)))^ 



+ ^ / \\Vh N (;u)\\ 2 du\F s 



N 



AV N (0). (32) 
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By assumption (|20|) . 

AFi(O) = A(Wi * Wi)(0) = (VWi * Wi)(0) 
VWi(-y)VVFi(y)dy 

< oo. 
As a consequence 

|AViv(aO| =N l3 N 2 l 3 ' d \AV l {N l3 l d x)\. 

So 2 

°^~*W (0) = Cl a%(t - s)N^'^\ 
and the thesis follows. 



(33) 

(34) 
(35) 



Remark 0.1. A^it) is not negative; indeed in general 

9 , ,9 a 2 , b 2 , a b . 
a 2 -ab + b 2 > ab + — = (— -) 2 > 0. 

Let now consider a special class of test functions, i.e. positive function 
6 G C 2 (M d ) such that 



4>{x) = \x\ for |a;| > 1 and 1 1 V0| U + | |A0| ^ < oo. (36) 
Lemma 0.2. For < s < t < T and (f> G C 2 (M d ) suc/i tfta* (E6j) ZioWs 

(.Xjv(i), 0) + C2-Aat(£) + zs a submartingale (37) 

and 



(Xj\r(i),0) — C2-Ajv(£) — is a supermartingale, 
with C2,C3 G R_|_ . 

Proo/. 

By applying Ito's Formula to (X/v(i),0), 
E[(Xjv(t),$|.F a ] 



(38) 



{x N (s),4>) 



r 



+E 



<Xjv(«), (-W(-) + (VG a * X JV (u))(0 - Vg N (;u))V(t> 



+^Acj))du\F s 
> (X N (s),<t>) 



-c 2 E 



-Vtf(-) + (VG a * Xjv(«))(-) - Vg N (-,u) + l)du|Jv, 



(39) 



Since 



< {X N (u),\-VU(-) + (VG a *X N (u))(-)-Vg N (;u)-l\ 2 ) 
= (X N (u), | - W(-) + (VG a * X N (u))(-) - Vg N (-,u)\ 2 
- 2(-VU(-) + (VG a * X N (u))(-) - Vg N (-,u)) + 1), 

2(X N (u), -VU(-) + (VG a * X N (u))(-) - Vg N (; u)) 

< (X N (u),\ - VU(-) + (VG a * X N (u))(-) -Vg N (-,u)\ 2 + 1), 

and therefore 

(X N (u), -VU(-) + (VG a * X N (u))(-) - Vg N (; u)) 

< (X N (u), | - W(-) + {VG a * X N (u))(-) - Vg N (-,u)\ 2 + 1). 

This implies that (|39p is greater than or equal to 

(X N (u), | - VU(-) + (VG a * X N (u)){-) - Vg N (;u)\ r < 



(X N (u), | - VU(-) + (VG a * X N (u))(-) - Vg N (-,u)\' 



(X N (s),<P) -c 2 E 
+2)du\F s ] 
> {X N (s),<j>) -c 3 E 
+l)du\F,] 

>(X N (s),<P)-c 3 E [ (X N (u),2(\-VU(-) + (VG a *X N (u))(-)\ 2 

-(-VU(-) + (VG a *X N (u))(-))Vg N (-,u) + \Vg N (-,u)\ 2 )) + a 2 N \\Vh N (-,u)\\ 2 du 

+ / du\!F s 

J s 

= {X N (s),<j>) - c 3 E [A N (t) - A N (s) + t - s\F s ] 

= (X N (s), <P) - c 3 E [A N (t) + t\F 8 ] + c 3 A N (s) + c 3 s. (40) 
Hence, 

E [(X N (t), <f>) + c 3 A N (t) + c 3 t\F s ] > (X N (s), <P) + c 3 A N (s) + c 3 s 



and ([37|) follows. 

In a completely analogous way (with +c 3 instead of —c 3 ), we obtain the 
property (|38l) . 



□ 



Let us define the sequence of stopped processes XN^(t) = Xj\f(tAT^), < 
t < T, N € N, and k > fixed, where rjy is defined by 

4 = inf{i>0 : HM-^lll + ^Mi) 

- / (X N (u),2\ - VU(-) + (VG a * X N (u))(-)\ 2 )du > k}, (41) 
Jo 

where k € M+. 

We consider a slight modification of the more general characterization of 
the relative compactness by Ethier-Kurtz (see [8], theorem 8.2). In particular 
we prove confining compactness property for the process {X/v,fc(i)} and then 
the boundedness of small variations of the process. 
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Proposition 0.1. For any e > there exists a compact K k in (Ai-p(M. d ), \\-\\i] 
such that 

inf ¥{X Nk (t) G K k , Vt G [0,21} > 1 - e. 

Proof. 

Let Bj = {x £ M d : \x\ > A, A > 1}; 
(X Ntk (t),<f>) = [ <f>(x)X N , k (t)(dx) > 



{x:\x\>\, A>1} 



<j>(x)X Ntk (t)(dx) 



\x\X N>k (t)(dx) > A / X Ntk (t)(dx) 

{x:\x\>\, A>1} ./{x:|x|>A, A>1} 

= A(Xjv )fc (i), Isj); 
if (Xjv,fc(t), 1b c ) > 5, then (X^^{t)-,4>) > A<5 and therefore 



sup(Xjv fc (t), Ibf) > <5 



sup(X7v,fc(i),0) > A<5 



< P jsup(A> )fc (t), 4>) + c 3 A N (t A 4) + (t A t£)c 3 > A<5 j . (42) 

By (|37|) in Lemma 1(01 and by Doob's Inequality, ([4"2"]1 is less than or equal to 
1 



A<5 



V[(X Njk (T),<l>)]+C3E A N (TArfr) 



+ c 3 E 



1 



< — [E[E[(X Ntk (T),(j>)\F ]] + c 3 E A N (TAr k Nj 



T A t% 



+ c 3 E 



T At 



N 



a i 

< — E 
" A<5 



{X Ntk (0), 4>) + c 3 E[A N (T A t^)\T ] + c 3 E[T A T k N \T Q ) 
c 3 ^(0 A 4) - (0 A T k N )c 3 + c 3 Ajv(T A 4) + (T A r&)c 3 ] ) . (43) 



By the definition of and since VG a , V£7 G 



Atv(tat^) < fc-IIM-.tA^nl 



TA- 



- / {X N (u),2\ - W(-) + (VG a * X,v(«))(-)n^ 
'o 

< k + / (X^(n),2| - VC/(-) + (VG a *X^(n))(-)| 2 )^ 

J o 

< + c 4 T; (44) 



it follows that (J43J) is less than or equal to 

1 JUJ c *(k T) 

— (E[(X N ^ k (0)^)+c 3 (k + c 4 T) + c 3 T + c 3 (k + c 4 T) + c 3 T]) < J±-L-l.(45) 

AO AO 



As a consequence, 



sup(X^(t),l B c) >( jl<^Ml. 



(46) 



in 



Let us now take e > and two sequences /ij and <5j of positive numbers such 

•■s(k,T 



that A*i = 6 and <*» \ °- Let X i = ^Bp- ~> 00 • Then © yelds 



Pjsup(Xjv, fc (*),%.) ViGNj < ^Pjsu^X^^),^.) >Sij 

00 n m\ 00 
(=1 1=1 

By Prohorov's Theorem, the set 

K k t = {fx E AMlR d ) = (M. 1b? Vi E N} 

is compact in A^-p(M rf ); since 

pjsup^fc^.lfl^) ><5 4 ,Vi E n| = l-p{(Xtf jfc (i),lBc.) < Ji,V» E N,V< E [0,T]} , 

by (|4"7)). Ve > there exists a compact set K k C .Mp(M d ) such that 
inf ¥{X Nk (t) E K k , Mt E [0,T]} > 1 - e. 

N&i 

Next proposition states that for little time variations we have little varia- 
tions of the processes {Xjy^it)}. 

Proposition 0.2. For any < 5 < 1, there exists a sequence {7^(<5)} ne N of 
non negative random variables such that 

E[\\X Nik (t + S)-X Njk (t)\\f]<E[^{S)] 0<t<T (48) 

and 



□ 



lim limsupE[7^(<5)] = 0. (49) 



Proof. 



1 I 



N 



\\X N , k (t) - XxMWi = sup - £ {/(Xjf^t)) - /pfyfcto)) 

1 * 

<^Ei^(*)- x ^( s )i 

i=l 

1 * 



i=l 



AT 

i=l 



-VE/p^(u)) + (VG a * ^(«))(^(«)) - Vg N (X i N (u),u)du 



+ <J N (W l {tf\T k N )-W i {sf\T k N )) 



C7 



N 



<-Y 

N ^ 

i=i 

N 



tAri. 



(ill 



+ 



N ^ 



N 



|(VG a * Xjv(«))(^(u)) - V 5A r(X^(n),n)| 



i=l 



(50) 



£7 



+jEK l ("^)-nA4) 



£ 

i=l 

A' 



f tAr N 






/ 


+ 




JsAt^ 




IsAt% 



(51) 



8=1 



By the Cauchy-Schwartz and Jensen inequalities, 



sAr; 



(Xjv(n), |(VG a * X N {u)){-) - Vg N (-,u)\)du 



< 



< 



tArfc \ X l 2 ( rtAr 



du 



sAt n J 



(X N (u), |(VG a * X N (u))(-) - Vg N (-,u)\) 2 du 



\ 1/2 



<Ar? 



\ 1/2 



du 



sAtJ 



<AtJ 



(X^n), |(VG a * - Vg N (; u)\ 2 )du 



1/2 



(52) 



moreover, if s < r^r, 



and if s > r^r 



tAr<? 



sAt! 



(in 



(At 



tin < / du = t — s 



tArJ 



(in 



sAt 



(in = < t — s. 



(53) 



(54) 
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Therefore, by (|52p ,(|53 |) and ((Ml), dSQJ) is less than or equal to 

*X N (u))(-) - Vg N (-,u)\ 2 )du 



( f tAT N 

c 7 \t - s\ + (t - s) 1 / 2 I J (X N (u),\(VG Q 



+^EK^ Ar £)-^( sAr ^ 



N 



N 



i=l 



( r tAT N 

< c 7 \t - s\ + (t - s) 1 / 2 ( J (X N (u),2\(VG a *X N (u))(-)f 



1/2 



-2((VG a *X N (u))(-)Vg N (',u)) + 2\Vg N (',u)\' 

+a 2 N \\Vh N (-,TAT k N )\\ 2 du+\\h N (-,TAr^)\\ 2 
N 

1=1 



c 7 \t -s\ + (t- s) 1 ' 2 (A N (T A r k N ) + \\h N (; T A r&)| ||) 



1/2 



A r 



i=l 



< c 7 \t - s\ + (t - s) 1 ' 2 ik + J (X N (s),2\(VG a 



1/2 



* X N (u)){-)\ 2 )du 



+^EK^ A ^)-^( sA ^) 



JV 



iV 



i=l 



< c 7 \t -s\ + (t- sfl 2 {k + c 6 T) 1 ' 2 + ^Y1 \ W ^ A t n) ~ W '( s A t n) 



N 



i=l 



(55) 



It follows that 
\\X Ntk (t)-X Ntk (s)\\i<2 

As a consequence, for < s < t < T 
E[\\X Njk (t)-X Ntk (8)\\i] < 2 6 E 

in particular, with t — s = 5 and 
7^ = S ^ 2 6 

we obtain 



4(t - s) 4 + (t - s) 2 (fc + c 6 T) 2 + ^(i - s) 2 



(56) 



4(t - s) 4 + (t - s) 2 (k + c 6 T) 2 + ^L(t - s) 2 



c^ + 5 2 (k + c 6 T) 2 + ^5 2 



E [\\X Njk (t) - X Njk (s)\\f\ < E[ 7 £(«5)] 
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and 



lim supE[ 7 ^(<5)] = 0. 



Proposition 0.3. {C{X^{- A Tn))}n£N, the sequence of probability laws 
of the processes {X^((t A rjy)), < t < T} is relatively compact in 
M r (C([0,T],Mv(R d ))). 

Proof. 

It is an obvious consequence of Proposition 10.21 and Theorem 8.6 p. 137, in [8]. 



Proposition 0.4. For any r such that < r < oo, 



lim inf P{r£v > rj = 1. 



Proof. 

By Lemma [Q.lj, the process 



t^ S N (t) = \\h N (;t)\\ 2 +A N (t)- f \x N (u),2\-VU(-)+(VG a *X N ( U ))(-)\ 2 )du 

Jo 

is a submartingale. 

By Doob's inequality 

sup S N (t) >k\< ~E[SW(t)] = \nM N {r) + Ta 2 N N^ d+2 ^ d - l Cl ] 



t<T 



E[M n (t)\ Fo] + raj.N^y^d 
yE[M N {0) +ra 2 N N^ d+2 ^ d ~ 1 c 1 ] 



\[n\hN{-M\l]+ra 2 N N^ 2 )/^ Cl 



(57) 



since limjv-»oo °~N = °~oo > 0, by (|2"2j) and ([25]) or ([501) . ([57]) is less than or 
equal to cg(r)/k, uniformly in N. 
The thesis follows. 



Remark 0.2. Proposition \0.4\ implies that t A rjy = t, for any r such that 

0<t<T. 

Proof of Theorem \U.1[ 

At this point Theorem 10.11 simply follows from Propositions 10.31 10.41 and Re- 
mark F 



1 J 



Theorem 10.11 implies the existence of a subsequence Nk C N, 
Ni < N2 < such that the sequence {£(X/v fc )}fceN converges in 

M-p(C([0, T],M-p(M. d ))) to some limit C(X), which is the distribution of some 
process {X(t),t G [0,T]}, with trajectories in C([0,T\, Mp(R d )). We discuss 
the uniqueness of the limit later on. By now we assume the uniqueness, so 
that {N k } = N. 

By Skorokhod's Theorem, we are allowed to assume that {Xpf(t),t G 
[0,T]} converges P-almost surely to {X(t),t G [0,T]} as N grows to infin- 
ity. So, we have 



lim sup\\X N (t) - X(t)\\i = 

N->oo t< T 



a.s. 



(58) 



Absolute continuity of the limit 



Next proposition deals with the regularity properties of the limit measure 
X(t). We consider the viscous case liniTv^oo &n = &00 > 0. 

Proposition 0.5. Suppose that limTv^ooOAr = > 0. For any t > 0, the 

measure X(t) is absolutely continuous with respect to Lebesgue measure on R d 
with a density p G L 2 (R d , R + ). 

Proof. 

We begin by showing that there exists a positive function p(x, t) such that 

-T 

0. 



lim E 

N-*oo 



lim E 



lim E 

N,N'-^oo 



< lim E 

N,N'^oo 



JR d 
T 



JM d 
T 



JR d 
T 



\hn{x,t) — p(x,t)\ 2 dxdt 

\h]sr(x,t) — hN'(x,t)\ 2 dxdt 
\h^(X,t) - h^>{X,t)\ 2 dXdt 



(59) 



J{\X\<k} 
T 



\h N (X,t) - h N >(X,t)\ 2 dXdt 



+2 lim E 

N,N'->oo 



\h N (X,t)\ 2 + \h N ,(X,t)\ 2 dXdt ; (60) 

'0 J{\X\>k} 

since, by (gDJ), |W^(A)| is bounded and ^v(A,t) = (X N (t),e iX -)W^(X), ex- 
pression (|6"0]) is less than or equal to 

\{X N (t),c° 

l0 ^{|A|<fe} 

rT r |\|2 |\|2 

l -^-\h N (X,t)\ 2 + i-J-|/uv(A,t)| 2 dAdt 



lim E 

7V,iV'^oo 



\{X N (t),e lX -) - {X N ,(t),e lX -)\ 2 dXdt 



+2 lim E 

N,N'^oo 



J{|A|>fc} 



A- 



A- 



< lim E 



T 



\(X N (t),e iX -) - (X N ,(t) : e lA -)\' z dXdt 



iX-\\2. 



J{\X\<k} 
T 



H — lim E 

k N,N'^oo 



Jw. d 



\X\ 2 \h N (X,t)\ 2 + \X\ 2 \h N ,(X,t)\ 2 dXdt 



(61) 
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Since by {58 



lim sup sup \(X N (t),e iX -) - (X N ,(t), e iX ')\ = 0, VA: > 

7V,AT'->oo t < T \x\<k 



expression (j6T|) is equal to 



a.s., 



Now 



- lim E 

k N,N'^oo 



E 
E 



|A| 2 |/ijv(A,t)| 2 + \\\ 2 \h N >(\,t)\ 2 d\dt 



T [ \\\ 2 \h N (\,t)\ 2 + \\\ 2 \h N >{\,t)\ 2 d\dt 

JR d 
T r T 

2j, , / \\T7U--.f *M|2 



(62) 



\\Vh N (;t)\\ dt+ / \\Vh N ,(;t)\\dt 



(63) 



by ([57|) . with A]y(t) as denned in (|3T|) and S]y(t) denned in Proposition [0? 
we obtain that 



E[||VM-,t)H 2 ]<^# 1 + ^ <- 



'JV 



a 



(64) 



jv 



uniformly in N with c positive constant. As a consequence (|63h is finite. 

It follows that, for sufficiently large, (|62p can be made smaller than any 
given e > and there exists a positive function p(x, t) G L 2 (M d , ]R + ) satisfying 



equation (|59j) . 

Since by (f59l) and lirmv-xx) Wjv(-) = 5q (in the sense of distributions) 



lim / f(x)X N (t)(dx) 



f(x)p(x,t)dx /£C 6 Vx[0,T]) 



a.s., 



we have by 

f(x,t)X(t)(dx) 



f(x,t)p(x,t)dx f G C£(M d x [0,T]),P- a.s. 



Therefore the measure X(t) is absolutely continuous with respect to the 
Lebesgue measure with density p(x,t). 

As a consequence of Proposition (|0.5|) 

lim (X N (t), /(•)) = /(•))= / f(x)p(x,t)dx /£C°(4te[0,T] 

(65) 

As next point we need the description of the dynamics governing the time 
evolution of the possible limit process {X(t),t G [0,T]}. 



□ 



A formal derivation of the continuum mod- 
els 

In this section, following [12], we characterize the limit behavior, as N — > oo, of 
the process Xn both in the case liniAr_ too on = and lim/vr^oo on = <7oo > 0. 



1 fi 



By taking into account expression (|14p and by using Ito's formula we get 
the following weak form of the time evolution of Xjy(i): 

(X N (t),f(;t)) = (X N (0),f(;0)) 

+ f (X N (s),(VG a *X N (s))-Vf(;s))ds 
\x N (s),Vg N (;s)-Vf(; S ))ds 
\x N ( S ),VU(-)-Vf(-,s))ds 
j\x N (s), ^a%Af(; s) + ^/(., s))ds 

I* N 

1/0 fc=l 



+ 



x [0,T]). 
(66) 



Last term in ([66 



M N (f,t):=^ ^Vf(X k N (s),s)dW k (s) 
Jo k=i 

is a martingale with respect to the natural filtration of the process {X]\f(t),t G 
[0, T]} and the quadratic variation 



lim E 

AT— >oo 



sup I M N (f,t) I 

t<T 



2 

= (67) 



(see [12], |14j). This implies, in both cases, convergence to zero in probability, 
that is the substantial reason of the deterministic limiting behavior of the 
process, as TV — ► 00, since in this limit the evolution equation of the process 
will not contain the Brownian noise anymore (see |14|). 

In order to derive a formal limit for the process Xn also when 
liniTv^ooCTV = 0, let us assume that X(t) admits density with respect to 
the Lebesgue measure also in this case. As a formal consequence of this as- 
sumption and ([65]) . we get 

lim g N (x,t) = lim (V N * X N (t))(x) = p(x,t), 
lim V<7Ar(x,t) = Vp(x,t), 

./V— >oo 

lim (VG a *X N (t))(x) = (VG a *X(t))(x) 

= jvG a (x-y)p(y,t)dy, x G M. d , t G [0, T]. 
Hence by applying the above limits, from (|66p an d the hypothesis 
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liniTv^oo ctv = <7oo > 0, we get the following equation 



f(x,t)p(x,t)dx = / f(x,0)p(x,0)dx 
jR d 
t 







+ ds [(VG a * p(-, s))(x) - VU(x) - Vp(x, s)) 



•V/(x, s)p(x, s)dx 



rt 

+ I ds 



o Jm. d . 



■2 



dx. 

(68) 



Equation (|68l) is the weak version of the following equation for the spatial 
density p: 

d a 2 

Q- t P( x ^) = ~f^p{x,t) + V -(p{x,t)Vp(x,t))+V ■ (p(x,t)VU(x)) 

- V- [p(x,t)(VG a *p(-,t))(x)], x G M d , t G [0, T], 

p(a;,0) = poW, x G M. d . (69) 



Obviously if (Too = the diffusive term in (I68j) and (|69p vanishes, while if 
o"oo > the dynamics of the density is smoothed by the diffusive term. This is 
due to the memory of the fluctuations existing when the number of particles 
./V is finite. 

Main results 

In the present section we present the main results of this chapter, namely 
two theorems on the convergence of the interacting particle system to the 
integro-differential equation (f69l) . both for = and > 0. 

We begin with the = (non-viscous case) following the approach 

proposed in [12] and then we move to the case <Too > (viscous case). 

Non- viscous case 

We are not aware of general results concerning the existence of sufficiently 
regular solutions p for equation (|69h : therefore we need the following assump- 
tion: 

Assumption 0.1. For some T G [0, oo) system 169\) with o"oo = admits a 
unique, nonnegative solution p G C^ rf+2 ^ //rf ] +1 ' 1 (IR d x [0, T]). 

About the uniqueness of the solution of equation (|69p without confining 
potential we address to [4]. 

Let (Too = and suppose that 

VW l G L°°{R d ), (70) 



\W!(x)\< d+2 , if M>1. (71) 



IS 



Consider the following assumption for the aggregation kernel G a (x) and 
the confining potential U(x): 

VG a ,VU e c l b {d+2)/d]+1 (R d ) n (72) 

As far as (3 is concerned, we need to assume that one of the following 
conditions is satisfied: 

< /J < jk- < 73 > 

or 

d 



< /3 < 1 and lim a N N^ d+2) / d - 1 = 0. (74) 



d + 2 N^+oo 

Under previous hypotheses, we can prove the following theorem: 
Theorem 0.2. Assume £7U\)-fi7$) and Assumption \0.1[ If 

lim E[\\h N (;0)-p (-)\\ 2 2 ] =0, 

iv— *-oo 



then 

0, (75) 



sup\\h N (-,t) - p(;t)\\l 
t<T 



lim E 

N^oo 

where p is the unique solution of WEfy with <7oo = 0. 
Corollary 0.1. Equation ( f75[ ) implies 

lim (X N (t),f) = (X(t),f) = f f(x)p(x,t)dx 

uniformly in t G [0,T], /or any f £ C£(R d ) n L 2 (M d ). 
Proo/. 

|(^(t)-p(-,t),/>| < K/iAr(-,t) - t), /)| + <JCiv(t), 1/ - / * Wiv|>; 
from (I75p and Lemma 10.31 we obtains our thesis. 



Previous corollary state that the empirical measure X]^(t) converges 
weakly to the density p(-,t). 

Proof of Theorem 10.21 

We prove this result by following the same approach used in [12] for proving 
it in the case of system ([9]) without confining potential (71 = 0). 

We obtain the convergence of /i7v(-,i) to its limit by performing the fol- 
lowing steps: 



1. we guess the dynamics (I69p for the limit p\ 

2. we try to control ||/i/v('> p('-> ^)l II m term of its initial value | |/ijv(-, 0) — 
Po(")lll by writing down Ito's formula for that process and by estimating 
the different contributions. 



□ 



1Q 



We have 

|2 lit / .\||2 , ii / .\||2 



\\h N (;t) - p(.,t)\\i = \\h N (;t)\\i + \\ P (;t)\\i - 2(p(;t),h N (;t)). (76) 

From (|68p and integration by parts, one gets: 



\\p(;t)\\l = \\P(-M\ 2 



2 

/ ds I [{VG a *p(-,s)){x)-Vp(x,s)-VU(x)]-Vp(x,s)p(x,i 
Jm d 



o 

IIp(-,o)|| 2 



ds 

2 



t 

ds I V • [p(x,s)(VG a * p(-,s))(x)]p(x,s)dx 



o 



+ ds V • [V p(x, s)p(x, s)]p(x, s)dx 
Jo JR d 

+ ds V • [VU(x)p(x,s)]p(x,s)dx 
Jo Jtt d 

+ ds V ■[Vp{x,s)p{x,s)]p(x,s) +V -[VU{x)p(x,s)]p{x,s) 
Jo JM. d 

- V • [p(x, s)(VG a * p(-, s))(x)]p(x, s)dx 

= \\p(-M\1 

— 2 ds V • [p(x, s)(VG a * p(-, s))(x)]p(x, s)dx 
Jo JR d 

+ 2 ds V • [V p(x, s)p(x, s)]p(x, s)dx 
Jo Jm d 

+ 2 ds V • [VU{x)p(x,s)]p{x,s)dx 

Jo JR d 



||p(-,0)||i+2 / (p(;s),(p(;s)*VG a )-V P (;s)}ds 





t ft 

2\ 



2 (p(;s),\V P (;s)\ 2 )ds-2 / (p(;s),VU(-)-Vp(; S ))ds. (77) 

Jo Jo 



With the same computations used to obtain expression (|32p . for the first term 



9fl 



of (|76p one obtains 

\\h N (-,t)\\ 2 2 = ||M-,0)||!-2 [\x N (u),\Vg N (; S )\ 2 )du 

Jo 

- 2 [ {X N (u),Vg N (;u)-VU)du 
Jo 
ft 



+ 2/ (X N (u),Vg N {-,u)-{VG a *X N (u)))du 
Jo 



on I \\Vh N (;u)\\idu 

t N 







+ / J2Vg N (X k N (s),s)dW k (s) 

Jo k=i 



+ c 1 a 2 N (t-s)N^ d+2 y d - 1 . (78) 

From (|66p . ()69|) . the symmetry of Wn, by Ito's formula and integration by 
parts, one obtains 

( P (;t),h N (;t)} = (X N (t), P (;t)*W N ) 

= <p(-,o),M-,o)> 

+ / <Xjv(t),[(Xjv(a)*VG a )-V^a)-V^(a:)]-V(p(-,t)*Wjv))d« 







2 rt rt a 

+ ^ J (X N (s),Ap(; 8 )*W N )ds + J^(X N (s),—p(; S )*W N )ds 



= (p(-,o),m-»o)) 

+ / (^(^[(XjvW + VGoJ-V^C-.flJ-V^x)]- Vp(-,t)*^jv)ds 

JO 

+-f / (Xtf(*),Ap(-,s)*W^>d« 

+ / (X 7V ( S ),[V-(p(-, S )Vp(-, S )) + V-(p(-, S )V[/) 

JO 

-VX-, fl )(VG a *p(-, a ))]]*Wiv)ds + -^ / ^V(j>(;s)*W N )(Xk(s))dW k (s) 

^° fc=i 
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(p(-,o),m-,o)) 



+ / (X N (t), [{X N (s) * VG a ) - Vg N {; s) - VU(x)] • Vp(-, t) * W N )ds 
Jo 

[\vh N (;s),Vp(;s))ds 

1 Jo 

+[ (Vh N (; S), [-Vp(-, S)-VU+ (VG a * P (; S))] P (; S ))ds 

Jo 

ft N 

+JF ^V(p(-, S )*^)(X^( S ))^ fc ( S ). (79) 
Jo k=i 

It follows that for (1T6|) one gets the following expression: 



\\h N (;t)-p(;t)\\ 2 



2 

HM-,o)-p(-,o)||l 

2 / [(X N (s), V 9N (; s) ■ (Vg N (; s) - Vp(; s) * W N )) 
Jo 



+ (p{; s),Vp(; s) ■ [Vp(; s) - Vh N (; a)])] ds 

+2 [ [(X N (s),((X N (s) *VG a ) -VU)[Vg N (;s) -Vp(-,s) *W N }) 
Jo 

+ (p(., a), ((p(; S) * VG a ) - VU)[Vp(; s) - Vh N (; s)})] ds 

-a% f [\\Vh N (;s)\\ 2 2 - (Vh N (;s),Vp(;s))] ds + C^t - s ) N Kd+2)/d-X 

Jo 

+2 ir I Y,( V 9N(;s)-Vp(;s)*W N )(X k N (s))dW k (s). (80) 



fc=l 



Estimates for the terms on the right side of ( 1801) 

In order to get estimates for the terms on the right hand side of (|80f) we need 
the following lemmas: 



Lemma 0.3. Let W\ be a symmetric density function which satisfies {70\j 
and (7ZJ) and let W N be defined as in {77p. Let f, V/ E C£(R d ). Then for 



any i£l 

\f{x)-{f*W N ){x)\<C 2X - N 1 \\yf\\oo- 

Furthermore if /, V/ G L 2 (R d ), then 

\\f-f*W N \\ 2 2 <c 2X ], 2 \\Vf\\l 
The constant c 2 is independent of f. 

Lemma 10.31 is proved in [17] . 
Lemma 0.4. Let us suppose Assumption \0. 1\ and let be L = [{d + 2)/2] and 

v(x) G C^ +1 (R d ,R d ), 
G{x) G L 1 (R (i ,R). 
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Then we have 

\(X N ( S ) - p(; 8), (V[h N (; S ) - p(; S )] * W N ) ■ v) \ 

<c 3 (\\h N ( S )-p(., S )\\l + X - N " + N^-^), 

| ((X N ( S ) - p(; S )) * G, (V[h N (; s) - P (; s)} * W N ) • v)\ 

<c 4 (\\h N ( S )-p(.,s)\\l+ X ~ N 2 + N^-^). 
For a proof of this lemma see [12] . 

We begin by considering the second term on the right hand side of (|8U|) ; 
since this term does not depend on the potential U, we can recall a result 
proved in [12]. Indeed, as showed in [12], by considering v = Vp(-,s) * Wn, 
from Lemma 10. 3| Lemma 10.41 and Assumption 10. 1[ we have 



A := \{X N (s),Vg N (;s)-(Vg N (;s)-Vp(;s)*W N )) 

+ (p(;s),Vp(;s)-[V P (;s)-Vh N (;s)})\ 
< C 5 \\h N (;s)-p(;s)\\ 2 2 + XN 2 + N P{1 " 2L/d) . (81) 

Now let us consider the term in (|80p involving the aggregating kernel G a 
and the confining potential U 

B = {X N (s),((X N (s)*VG a )-VU)-[V(h N (;s)-p(;s))*W N }) 

- (p(; s), ((p(; s) * VG a ) - VU) • [~V p(- , s) + Vh N (; S )]) 

= {X N ( S ) - p(; S), ((X N (S) * VG a ) - W) ■ [V(h N (; s) - p(; S )) * W N ]) 

+ (p(; s), ((X N (s) * VG a ) - VU) ■ [V(h N (-, s) - p(; s)) * W N ] 

-((p(-, s) * VG a ) - VU) ■ [-Vp(; s) + Vh N (; s)}) 
= B 1 + B 2 . (82) 

By ([72]) and Lemma E3 with v = -VU + X N (s) * VG a we get 

< c 3 (||M*) - p(; s)\\l + x~ N 2 + N^- 2L ^) . (83) 
On the other hand 

B 2 = (p(-, s), {(X N (s) * VG a ) ~ VU) • [V(h N (-, s) - p(-, s)) * W N ] 
-((/>(■, s) * VG a ) - VU) ■ [-Vp(; s) + Vh N (; a)]) 

= (p(; S), [(X N (S) - p(; S )) * VG a ] • [V(h N (; s) - p(; S )) * W N }) 
+ (p(;s),((p(;s)*VGa)-VU) 

V[(h N (; s) - P (; «)) * W N - (h N (; s) - P (; s))}) 
= ((X N ( S ) - p(; S)) * VG a ,p(; s)[V(h N (; s) - p(; S )) * W N ]) 
-(V[p(; S )((p(;s)*VG a )-VU)], 

(h N (-, s) - p(-, s)) *W N - (h N (-, s) - p(-, s))) 
= B\ + Bl (84) 
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By Assumption 10.11 condition (|72p and Lemma 10.41 with G = VG a and 
v = p — VU we get 

\B\\ < C4 (||M", 8) ~ p(; S)\\l + X J* + NfiW>) 

and by taking into account also Lemma 10.31 

|B 2 2 | <c 2 (||M^)-pM||| + x^ 2 ). 

So for the second term B of (|80jl we get the following estimate: 

\B\ < c 6 (| |^(, s) - p(; s)\\l + x~ N 2 + NM-W) . 
For the third integrand on the right side of (|80j) we have 

C := -a 2 N [\\Vh N (;s)\\ 2 2 -(Vh N (;s),Vp(;s))} 
'\\VhsM\\l 



< 



2 



(Vh N (;s),Vp(;s)) 



-^L\\Vh N (;s) - Vp(;s)\\ 2 2 + ^L\\Vp(; S )\\l 



Let us consider the sub martingale term in (|80j) : 

rt N 
'° k=l 

As showed in [12], by Doob's inequality, 

rT' 



(85) 



N 



J2(V9n(; s) - Vp(-, s) * Wtf)(X&(s))dW fc (a) 



E 



sup Mjv(t) 

t<T' 



< c 7 a 2 N x d N E 
N 



\\V(h N (.,s)-p(;s))\\ 2 2 ds 







T' < T. 

(86) 



By collecting all contributions, (|80p becomes 



||M-.<)-p(-.*)ll2 < ll^(-,0)-p(-,0)||^+C5 / ||/lAf(-,s) -p(-,s)||^S 

'0 



+ c 5 t[ X - N 2 + N^- 2L ^ + a 2 N N^^ + a 2 A 



a 



N 



\V(h N (-,s)- p(-,s))\\ 2 ds 



+ 



2<7jV 

N 



N 



Y,(V9n(; s) - Vp(; s) * W N )(X k N (s))dW k (s) 



k=l 



(87) 



for < t < T. 
From (IHBT) and 



E 



sup||M^)-p(^)Hl + ^(l-^ 

t<T> 1 \ iv 



\V(h N (;t) - p(;t))\\ 2 2 dt 



< E [11^0,0) -p(;0)\\l 



T' 



\\h N (;t) - p(.,t)\\ 2 dt 



o 

i3(d+2) 
d 
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For N sufficiently large and by applying Gronwall's inequality we obtain 



E 



sup||/iAr(-,i) - p{-,t)\\l 



t<T 



< 



E [\\h N (; 0) - p(, 0)||1] + T( Xn 2 + N PW> + *%N^-i + a%) 



As N — > oo, by (|74|) and since liniTv-^oo °at-^ d = 0, we obtain our thesis. 



Viscous case 

Now we move to the case o"oo > 0. 

Due to technical difficulties (the presence of the non vanishing term > 
), in this case we can not carry out the same proof as for Theorem 10.21 
Therefore, by following [15], we try to control directly E [(X(t), f) — (p(t), /)], 
obtaining a result analogous to Corollary 10.11 

About the regularity and uniqueness of the solution of equation (I69j) we 
make the following assumption 



Assumption 0.2. System I169\) with Goo > admits a unique, nonnegative 
solution p E C 2 > 1 (R d x [0,T]). 

The requirements of Assumption 10.21 are weaker than those of Assumption 
EE] (p E C 2 * 1 ^ x [0,T]) instead of p E C ^ d+2),d]+1 ' 1 (M. d x [0,T])), but we 
need a further restriction on the function W\ defined by (fT9j) : W\ must have 
compact support. 

Theorem 0.3. If 

i) 

lim C(X N (0)) = S M in M r (M P (R d )), (90) 

N— *oo 

where po has density p(x, 0) with respect the Lebesgue measure, 
ii) the parameter (3 satisfies condition \30}) . 
Hi) W\ defined in M9\j has compact support, 
then 

lim (X N (t)J(-,t)) = (X(t)J(-,t)} = f f(x,t)p(x,t)dx (91) 

for any f E C^' l (M. d ,M.+ ), where p is the unique solution of \69\) with > 0. 
Proof. 

We have to show that 

SUp \\{X(t),f{;t)) - I f(x,t)p( X ,t)dx\\ 1 = 0. 
0<t<T JW 



9£ 



□ 



Since (J68J) is the weak form of (|69p . it is sufficient to show that for any / £ 



E 



1 2 • . d 



(X(t), f(-,t)) - </i 0) f(; 0)) - J (p(; s), -o4A/(-, a) + ^/(-, a) 
+ [(VG a * *))(•) - W(-) - Vp(-, a)] • V/0, s))ds \] = 0. 



For fixed / G C 6 2,1 (M d , R+) 



E 



9 



{X(t),/(-,t))-{/i 0> /(-,0))-^{p(-, S ),^ c 



■-ai o Af(; S ) + -J(;S 



+ [(VG a * p{; a))Q - W(-) - Vp(-, a)] • V/(-, s))ds |] 



< E[|(X(*),/(.,i))-(X*(t),/(.,t))|] 
+ E[|( / uo,/(-,0))-(X JV (0),/(-,0))|] 



+ 



-E 



+ E 

+ E 

+ E 

+ E 



\(-p(; S ),Af(; S )) + (X N (s),Af(; S ))\ds 
(p(, S ),A /( , s)) + (XAr ( S ),^/(., S ))|d S 



| <p(-, s), Vp(-, a) • V/(-, a)} - (M"> a), V^(-, s) • V/(-, a)) |da 



<M", a), V%(-, a) • V/(-, s)) - {X N (s),Vg N ■ V/(-, s))ds 



-(p(; a), [(VG a * p(-, a))(-) - VCT(.)] • V/(-, a)) 



+ (X^a), [(VG a * X N (s))(-) - W(-)] • V/(-, s))\ds] 



+ E 



1/11 fe=i 



^( s ) )S )dWfe( s ; 



+ E[|(Xjv(t), /(•,*)) - (Xjv(0),/(-,0)) 

" (X Af ( s ),(VG *X JV (s))-V/(-,5))ds+ / (X JV ( s ),V^(-, S )-V/(-, S ))d S 

'0 JO 

+ J (X N (s),VU(-)-Vf(; S ))ds- J (X N (s),^a 2 N Af(;s) + -^-f(;s))ds 

rt N 



~N 



Y,Vf(X k N (s),s)dW k (s) 



k=l 



(92) 



i=l 



Clearly, by ([65]) and hypothesis ([90]), E^i^W = 0, by ([66]) I%(t) = 
and by (f67|) = 0. It remains to estimate the terms I%(t), I%(t) and 



9fi 



I 7 N (t). 



= E 

< IIA/Hoc / E 

J o 

< I|A/||oo(e 
r r T 

E 



\(p(; S ),p(;s)Af(;s)) - (h N (■ , s) , h N (■ , s) Af (■ , S ))\ds 



\h,N(x,t) - p(x,t)\\JiN(x,t) + p(x,t)\dx 



tit 



cT 




|/iat(x, t) — p(x, t)\ 2 dxdt 



1/2 



\Iin(x, t) + p(x, t)\ 2 dxdt 



10 



1/2 



by (|57|) and (|59|) we obtain 



lim I%(t) = 0. 

AT— >oo 



(93) 



By the symmetry of W\, 



E 



{X N (s),W N *(Vh N (-,s)-Vf(-,s)) 



iW N *Vh N (;s))-Vf(;s))ds\] 



E 



X N (s)(dx) / Wiv(s-i/)Vfeiv(y,s) 



<Vf(y)-Vf(x))dy) ds\] 



(94) 



By the definition of Wn and since Wi has compact support, with c 

i,j<d 



diam(suppWi(-)) and 1 1/? 2 /| |oo = supj j <d \ \d 2 A |oo, (fM)) is less than or equal 



to 

cX^II^/llooE 

< C x^||^ 2 /IIoc(e 

< CX^||^ 2 /Hoo. 

It follows that 



(X N ( S )*W N ,\Vh N (; S )\)ds 



T 



||^(-,s)||lds 



1/2 



E 



\\Vh N (-,s)\\ 2 2 ds 



1/2 



lim = 0. 

N— >oo 



(95) 



E 



-(p(-, S ),[(VG a *p(-, s ))(-)-VC/(-)]-V/(-, S )) 



+ (Xjv( S ), [(VG * Xat( s ))(.) - VC/(-)] • V/(-, s)) 
+ <X iV ( S ), [(VG a * p(; *))(•) - W(-)] • V/0, *)> 
- [(VG a * «))(■) - VU(-)] ■ Vf(; S))\ds] 



< E 



I (X N (s) - p(; s), [(yG a * p(; *))(•) - W(-)] • V/(-, a)) I 



+ \(X N (s), [(VG a * P(; *))(•) - (VG a * Xjv( s ))(.)] • V/(-, ds] 



97 



By Assumption [02] and ([58]) 



lim lUt) = 0. 



N— >oo 



As a consequence 




9 

£ 

i=l 



□ 



Conclusions 

In this paper we have studied the asymptotic behavior of system ([9]) for the 
size of the population N growing to infinity, being the time t fixed, in terms 
of a law of large numbers for the empirical process {X^(t),t £ 

It is also of interest to study the limiting behavior of such a system for 
fixed N and time growing to infinity. In [B] the authors investigate conditions 
for the existence of an invariant measure for system ([9]), i.e. conditions about 
the interaction potential and the confining potential such that there exists 
an invariant measure for the particle positions and, as a consequence, for the 
empirical process {X]\r(t),t £ 
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